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Abstract 

The Maximal Independent Set (MIS) problem is one of the basics in the study of locality in 
distributed graph algorithms. This paper presents an extremely simple randomized algorithm 
providing a near-optimal local complexity for this problem, which incidentally, when combined 
with some known techniques, also leads to a near-optimal global complexity. 

Classical MIS algorithms of Luby [STOC’85] and Alon, Babai and Itai [JALG’86] provide the 
global complexity guarantee that, with high probability^, all nodes terminate after 0(log n) rounds. 
In contrast, our initial focus is on the local complexity^ and our main contribution is to provide a 
very simple algorithm guaranteeing that each particular node v terminates after 0(logdeg(u) -I- 
logl/e) rounds, with probability at least 1 — e. The guarantee holds even if the randomness 
outside 2-hops neighborhood of v is determined adversarially. This degree-dependency is optimal, 
due to a lower bound of Kuhn, Moscibroda, and Wattenhofer [PODC’04]. 

Interestingly, this local complexity smoothly transitions to a global complexity: by adding 
techniques of Barenboim, Elkin, Pettie, and Schneider [FOCS’12; arXiv: 1202.1983v3], we^ get a 
randomized MIS algorithm with a high probability global complexity of 0(log A) -|- log")^ 

where A denotes the maximum degree. This improves over the 0(log^ A) -|- result of 

Barenboim et ah, and gets close to the n(min{log A, Vlog n}) lower bound of Kuhn et al. 

Corollaries include improved algorithms for MIS in graphs of upper-bounded arboricity, or 
lower-bounded girth, for Ruling Sets, for MIS in the Local Computation Algorithms (LCA) model, 
and a faster distributed algorithm for the Lovasz Local Lemma. 


^As standard, we use the phrase with high probability to indicate that an event has probability at least 1 — 1/n. 
^ quasi nanos, gigantium humeris insidentes 






1 Introduction and Related Work 


Locality sits at the heart of distributed computing theory and is studied in the medium of problems 
such as Maximal Independent Set (MIS), Maximal Matching (MM), and Coloring. Over time, MIS 
has been of special interest as the others reduce to it. The story can be traced back to the surveys of 
Valiant [Val82] and Cook [Coo83] in the early 80’s which mentioned MIS as an interesting problem 
in non-centralized computation, shortly after followed by (poly-)logarithmic algorithms of Karp and 
Wigderson [KW84], Luby [Lub85], and Alon, Babai, and Itai [ABI86]. Since then, this problem has 
been studied extensively. We refer the interested reader to [BEPSv3, Section 1.1], which provides a 
thorough and up to date review of the state of the art. 

In this article, we work with the standard distributed computation model called LOCAL [PelOO]: 
the network is abstracted as a graph G = {V,E) where |V| = n; initially each node only knows its 
neighbors; communications occur in synchronous rounds, where in each round nodes can exchange 
information only with their graph neighbors. 

In the LOCAL model, besides it’s practical application, the distributed computation time-bound 
has an intriguing purely graph-theoretic meaning: it identifies the radius up to which one needs 
to look to determine the output of each node, e.g., its color in a coloring. For instance, results 
of [Lub85, ABI86] imply that looking only at the 0(logn)-hop neighborhood suffices, w.h.p. 

1.1 Local Complexity 

Despite the local nature of the problem, classically the main focus has been on the global complexity, 
i.e., the time till all nodes terminate. Moreover, somewhat strikingly, the majority of the standard 
analysis also take a non-local approach: often one considers the whole graph and shows guarantees 
on how the algorithm makes a global progress towards it local objectives. A prominent example is 
the results of [Lub85, ABI86] where the analysis shows that per round, in expectation, half of the 
edges of the whole network get removed^, hence leading to the global complexity guarantee that after 
O(logn) rounds, with high probability, the algorithm terminates everywhere. See Section 2. 

This issue seemingly suggests a gap in our understanding of locality. The starting point in this 
paper is to question whether this global mentality is necessary for obtaining the tight bound^. That 
is, can we instead provide a tight bound using local analysis, i.e., an analysis that only looks at 
a node and some small neighborhood of it? To make the difference more sensible, let us imagine 
n ^ oo and seek time-guarantees independent of n. 

Of course this brings to mind locality-based lower bounds which at first glance can seem to imply 
a negative answer: Linial [Lin92] shows that even in a simple cycle graph, MIS needs D(log* n) 
rounds, and Kuhn, Moscibroda and Wattenhofer [KMWvl] prove that it requires Q{^/logn) rounds 
in some well-crafted graphs. But there is a catch: these lower bounds state that the time till all nodes 
terminate is at least so much. One can still ask, what if we want a time-guarantee for each single 
node instead of all nodes? While in the deterministic case these time-guarantees, called respectively 
local and global complexities, are equivalent, they can differ when the guarantee that is to be given is 
probabilistic, as is usual in randomized algorithms. Note that, the local complexity is quite a useful 
guarantee, even on its own. For instance, the fact that in a cycle, despite Linial’s beautiful D(log* n) 
lower bound, the vast majority of nodes are done within 0(1) rounds is a meaningful property and 
should not be ignored. To be concrete, our starting question now is: 

Local Complexity Question: How long does it take till each particular node v terminates, 
and knows whether it is in the (computed) MIS or not, with probability at least 1 — e? 


®These analysis do not provide any uniformity guarantee for the removed edges. 

^Without insisting on tightness, many straightforward (but weak) complexities can be given using local analysis. 
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Using A to denote the maximum degree, one can obtain answers such as 0(log^ A + log 1/e) 
rounds for Luby’s algorithm, or 0(log Aloglog A +logAlogl/e) rounds for the variant of Luby’s 
used by Barenboim, Elkin, Pettie, and Schneider [BEPSv3] and Chung, Pettie, and Su [CPS14]. 
However, both of these bounds seem to be off from the right answer; e.g., one cannot recover from 
these the standard O(logn) high probability global complexity bound. In the first bound, the first 
term is troublesome and in the latter, the second term becomes the bottleneck. In both, the high 
probability bound becomes 0(log^ n) when one sets A = for a constant 5 > 0. 

We present an extremely simple algorithm that overcomes this problem and provides a local 
complexity of 0(log A + log 1/e). More formally, we prove that: 

Theorem 1.1. There is a randomized distributed MIS algorithm for which, for each node v, the 
probability that v has not made its decision after the first 0(logdeg(u) + log 1/e) rounds is at most 
e. Furthermore, this holds even if the bits of randomness outside the 2-hops neighborhood of v are 
determined adversarially. 

The perhaps surprising fact that the bound only depends on the degree of node v, even allowing 
its neighbors to have infinite degree, demonstrates the truly local nature of this algorithm. The 
logarithmic degree-dependency in the bound is optimal, following a lower bound of Kuhn, Mosci- 
broda and Wattenhofer [KMWvl]: As indicated by [Kuhl5], with minor changes in the arguments 
of [KMWvl], one can prove that there are graphs in which, the time till each node v can know if it 
is in MIS or not with constant probability is at least H(logA) rounds. 

Finally, we note that the fact that the proof has a locality of 2-hops —meaning that the analysis 
only looks at the 2-hops neighborhood and particularly, that the guarantee relies only on the coin 
tosses within the 2-hops neighborhood of node v —will prove vital as we move to global complexity. 
This might be interesting for practical purposes as well. 

1.2 Global Complexity 

Notice that Theorem 1.1 easily recovers the standard result that after O(logn) rounds, w.h.p., all 
nodes have terminated, but now with a local analysis. In light of the H(min{log A, yTog^}) lower 
bound of Kuhn et al. [KMWvl], it is interesting to find the best possible upper bound, specially when 
log A = o(log n). The best known bound prior to this work was 0(log^ A) -|- rounds, due 

to Barenboim et al. [BEPSvS]. 

The overall plan is based on the following nice and natural intuition, which was used in the 
MIS results of Alon et al. [ARVX12] and Barenboim et al. [BEPSv3]. We note that this general 
strategy is often attributed to Beck, as he used it first in his breakthrough algorithmic version of 
the Lovasz Local Lemma [Bec91]. Applied to MIS, the intuition is that, when we run any of the 
usual randomized MIS algorithms, nodes get removed probabilistically more and more over time. If 
we run this base algorithm for a certain number of rounds, a graph shattering type of phenomena 
occurs. That is, after a certain time, what remains of the graph is a number of “small” components, 
where small might be in regard to size, (weak) diameter, the maximum size of some specially defined 
independent sets, or some other measure. Once the graph is shattered, one switches to a deterministic 
algorithm to finish off the problem in these remaining small components. 

Since we are considering graphs with max degree A, even ignoring the troubling probabilistic 
dependencies (which are actually rather important), a simplistic intuition based on Galton-Watson 
branching processes tells us that the graph shattering phenomena starts to show up around the 
time that the probability e of each node being left falls below 1/A®. Alon et al. [ARVX12] used 

®In truth, the probability threshold is 1/poly(A), because of some unavoidable dependencies. Bnt dne to the 
exponential concentration, the time to reach the l/poly(A) threshold is within a constant factor of that of the 1/A 
threshold. We will also need to establish some independence, which is not discussed here. See Section 4. 
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an argument of Parnas and Ron [PR07], showing that Luby’s algorithm reaches this threshold after 
0(Alog A) rounds. Barenboim et al. [BEPSvS] used a variant of Luby’s, with a small but clever 
modification, and showed that it reaches the threshold after O(log^A) rounds. As Barenboim et 
al. [BEPSvS] show, after the shattering, the remaining pieces can be solved deterministically, via the 
help of known deterministic MIS algorithms (and some other ideas), in log A • rounds. 

Thus, the overall complexity of [BEPSvS] is 0(log^ A)+log ^ = 0(log^ 

To improve this, instead of Luby’s, we use our new MIS algorithm as the base, which as Theo¬ 
rem 1.1 suggests, reaches the shattering threshold after only 0(log A) rounds. This will be formalized 
in Section 4. We will also use some minor modifications for the post-shattering phase to reduce it’s 
complexity from log A • to The overall result thus becomes: 

Theorem 1.2. There is a randomized distributed MIS algorithm that terminates after O(logA) -|- 
20 (Vioglogrt) with probability at least 1 — 1/n. 

This improves the best-known bound for MIS and gets close to the fI(min{log A, -^logn}) lower 
bound of Kuhn et al. [KMWvl], which at the very least, shows that the upper bound is provably 
optimal when log A E [2^^^°®^^^, yTogTi]. Besides that, the new result matches the lower bound 
in a stronger and much more instructive sense: as we will discuss in point (C2) below, it per¬ 
fectly pinpoints why the current lower bound techniques cannot prove a lower bound better than 
n(min{log A, y'log n}). 

1.3 Other Implications 

Despite its extreme simplicity, the new algorithm turns out to lead to several implications, when 
combined with some known results and/or techniques: 

(Cl) Combined with the finish-off phase results of Barenboim et al. [BEPSvS], we get MIS algorithms 
with complexity 0(log A) -|- 0(min{ -|- log A log log n, X + M log log n, A -|- (log log re)^"*"^}) for 
graphs with arboricity A. Moreover, combined with the low-arboricity to low-degree reduction 
of Barenboim et al. [BEPSvS], we get an MIS algorithm with complexity 0(logA -|- y/\ogn). 
These bounds improve over some results of [BEPSvS], Barenboim and Elkin [BEIO], and Lenzen 
and Wattenhofer [LWll]. 

(C2) The new results highlight the barrier of the current lower bound techniques. In the known 

locality-based lower bound arguments, including that of [KMWvl], to establish a T-round lower 

bound, it is necessary that within T rounds, each node sees only a tree. That is, each T-hops 

neighborhood must induce a tree, which implies that the girth must be at least 2r-|-l. Since any 

2 _ 

g'-girth graph has arboricity A < 0{n3-'^), from (Cl), we get an 0(vlogn)-round MIS algorithm 
when g = D(-y/log n). More precisely, for any graph with girth g = D(min{log A, -^logn}), we 
get an 0(min{log A-1-2*^^^^°®^ .^logn})-round algorithm. Hence, the D(min{log A, yTogTi}) 
lower bound of [KMWvl] is essentially the best-possible when the the topology seen by each 
node within the allowed time must be a tree. This means, to prove a better lower bound, 
one has to part with these “tree local-views” topologies. However, that gives rise to intricate 
challenges and actually, to the best of our knowledge, there is no distributed locality-based 
lower bound, in fact for any (local) problem, that does not rely on tree local-views. 

(C3) We get an 0{y/\og n)-round MIS algorithm for Erdos-R&yi random graphs G{n,p). This is 
because, if p = ^( 2'/*°^" ), with high probability the graph has diameter 0{y/logn) hops 

(see e.g. [CLOl]) and when p = 0( ^'^‘°^"" ), with high probability, A = and thus, the 

algorithm of Theorem 1.2 runs in at most 0{^^/logn) rounds. 
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(C4) Combined with a recursive sparsification method of Bisht et al. [BKP14], we get a (2, /3)-ruling- 
set algorithm with complexity 0(/3 log^'^^ A) + improving on the complexities 

of [BEPSvS] and [BKP14]. An (a,/3)-ruling set 5 is a set where each two nodes in S are at 
distance at least a, and each node v G V \ S has a node in S within its /3-hops. So, a (2,1)- 
ruling-set is simply an MIS. The term 0(/3 log^^^ A) is arguably (and even provably, as [Kuhl5] 
indicated,) best-possible for the current method, which roughly speaking works by computing 
the ruling set iteratively using /3 successive reductions of the degree. 

(C5) In the Local Computation Algorithms (LCA) model of Rubinfeld et al. [RTVXII] and Alon 
et al. [ARVXI2], we get improved bounds for computing MIS. Namely, the best-known time 
and space complexity improve from, respectively, log^ n and log^ n bounds 

of Levi, Rubinfeld and Yodpinyanee [LRY15] to log^n and log^n. 

(C6) We get a Weak-MIS algorithm with complexity O(logA), which thus improves the round 
complexity of the distributed algorithmic version of the Lovasz Local Lemma presented by 
Chung, Pettie, and Su [CPS141 from Oflog i re-log^ A) to Oflog i n-logA). Roughly 

ep(A-f 1) ep(A+l) 

speaking, a Weak-MIS computation should produce an independent set S such that for each 
node V, with probability at least 1 — 1/ poly (A), v is either in S or has a neighbor in S. 

(C7) We get an 0(log A -|- log log log n)-round MIS algorithm for the CONGESTED-CLIQUE model 
where per round, each node can send 0(logre)-bits to each of the other nodes (even those 
non-adjacent to it): After running the MIS algorithm of Theorem 1.1 for O(logA) rounds, 
w.h.p., if A > n^'^, we are already done, and otherwise, as Theorem 4.2 shows, all leftover 
components have size o{n^'^). In the latter case, using the algorithm of [HPP~*'I5], we can make 
all nodes know the leader of their component in O(logloglogn) rounds, and using Lenzen’s 
routing [Lenl3], we can make each leader learn the topology of its whole component, solve the 
related MIS problem locally, and send back the answers, all in 0(1) rounds. 


2 Warm Up: Local Analysis of Luby’s Algorithm 

As a warm up for the MIS algorithm of the next section, here, we briefly review Luby’s algorithm 
and present some local analysis for it. The main purpose is to point out the challenge in (tightly) 
analyzing the local complexity of Luby’s, which the algorithm of the next section tries to bypass. 

Luby’s Algorithm: The algorithm of [Lub85, ABI86] is as simple and clean as this: 

“/n each round, each node picks a random number^ uniformly from [0,1]; strict local 
minimas join the MIS, and get removed from the graph along with their neighbors! 

Note that each round of the algorithm can be easily implemented in 2 communication rounds on G, 
one for exchanging the random numbers and the other for informing neighbors of newly joined MIS 
nodes. Ignoring this 2 factor, in the sequel, by round we mean one round of the algorithm. 

Global Analysis: The standard method for analyzing Luby’s algorithm goes via looking at the 
whole graph, i.e., using a global view. See [MRIO, Section 12.3], [PelOO, Section 8.4], [Lyn96, Section 
4.5] for textbook treatments. We note that this is the only known way for proving that this algorithm 
terminates everywhere in O(logn) rounds with high probability. The base of the analysis is to show 
that per iteration, in expectation, at least half of the edges (of the whole remaining graph) get 
removed. Although the initial arguments in [Lub85,ABI86] were more lengthy, Yves et al. [YRSDZIO] 

®One can easily see that a precision of 0(log A) bits suffices. 
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pointed out a much simpler argument for this. See Appendix A, which describes (a paraphrased 
version of) their argument. By Markov’s inequality, this per-round halving implies that after 0(log n) 
rounds, the algorithm terminates everywhere, with high probability. 


2.1 Local Analysis: Take 1 

To analyze the algorithm in a local way, and to bound its local complexity, the natural idea is to say 
that over time, each local neighborhood gets “simplified”. Particularly, the hrst-order realization of 
this intuition would be to look at the degrees and argue that they shrink with time. The following 
standard observation is the base tool in this argument: 


Claim 2.1. Consider a node u at a particular round, let d{u) he its degree and dmax be the maximum 
degree among the nodes in the inclusive neighborhood N^{u) ofu. The probability that u is removed 
in this round is at least —. 

d[U)-\-dmax 


Proof. Let u* be the node in N~^{u) that draws the smallest random number. If u* actually has 
the smallest in its own neighborhood, then it will join MIS which means u gets removed. Since 
all numbers are iid random variables, and as u* is the smallest number of d{u) + 1 of them, the 
probability that it is the smallest both in its own neighborhood and the neighborhood of u is at least 
—. This is because, the latter is a set of size at most d(u) + dmax- □ 


From the claim, we get that if the degree of a node u is at least half of that of the max of its 
neighbors, then in one round, with probability at least 1/3, u gets removed. Thus, in a = 0(1) 
rounds from the start, either u is removed or its degree falls below A/2, with probability at least 
1/2. We would like to continue this argument and say that every 0(1) rounds, u’s degree shrinks 
by another 2 factor, thus getting a bound of O(logA). However, this is not straightforward as n’s 
degree drops might get delayed because of delays in the degree drops of u’s neighbors. The issue 
seems rather severe as the degree drops of different nodes can be positively correlated. 

Next, we explain a simple argument giving a weak but still local complexity of 0(log^'^A + 
logAlogl/e) rounds: For the purpose of this paragraph, let us say a removed node has degree 
0. From above, we get that after 10a log^'^ A rounds, the probability that u still has degree at 
least A/2 is at most Thus, using a union bound, we can say that with probability at 

least 1 — (A + 1)2“^^after 10a log^'^ A rounds, u and all its neighbors have degree at most 
A/2. Hence, with probability at least 1 — (A + 2)2“^°*°®^ '^^, after 20a log^'^ A rounds, node u has 
another drop and its degree is at most A/4. Continuing this argument pattern recursively for log®'® A 
iterations, we get that with probability at least 1 — (A + 2)*°s°'®^ • 2“^®^°®^ '^^ > 1 — 2“®*°®^ '^^, 
after lOalog^ A rounds, node u’s degree has dropped to A/2*°s Now, we can repeat a similar 
argument, but in blocks of 10a log^ A rounds, and each time expecting a degree drop of 2*°® ^ factor. 

We will be able to afford to continue this for log®'® A iterations and say that, after 10a log^'® A rounds, 
with probability at least 1 — (A + 2)*°®° * ^ • 2“®*°®^ * ^ > 1 — 2“^°®^ * the degree of u has dropped 
to 1/2. Since a degree less than 1/2 means degree 0, which in turn implies that v is removed, we 
get that V is removed after at most 0(log^'® A) rounds with probability at least 1 — 2“^^^°® A 
simple repetition argument proves that this generalizes to show that after 0(log^'® A + log A log 1/e) 
rounds, node u is removed with probability at least 1 — e. 

In the full version of this paper, we will present a stronger (but also much more complex) argument 
which proves a local complexity of 0(log^ A + log 1/e) for the same algorithm. This bound has the 
desirable additive log 1/e dependency on e but it is still far from the best possible bound, due to the 
hrst term. 
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2.2 Local Analysis: Take 2 

Here, we briefly explain the modification of Luby’s algorithm that Barenboim et al. [BEPSvS] use. 
The key is the following clever idea: they manually circumvent the problem of nodes having a lag in 
their degree drops, that is, they kick out nodes that their degree drops is lagging significantly out of 
the algorithm, as these nodes can create trouble for other nodes in their vicinity. 

Formally, they divide time into phases of 0(loglog A + logl/e) rounds and require that by the 
end of phase k, each node has degree at most A/2^. At the end of each phase, each node that has 
a degree higher than the allowed threshold is kicked out. The algorithm is run for log A phases. 
From Theorem 2.1, we can see that the probability that a node that has survived up to phase i — 1 
gets kicked out in phase i is at most Hence, the probability that a 

given node v gets kicked out in one of the log A phases is at most e. This means, by the end of 
0(log A log log A + log A log 1/e) rounds, with probability 1 — e, node v is not kicked out and is thus 
removed because of having degree 0. That is, it joined or has a neighbor in the MIS. 

This 0(log A log log A + log A log 1/e) local complexity has an improved A-dependency (and the 
guarantee has some nice independence type of properties). However, as mentioned in Section 1.1, its 
e-dependency is not desirable, due to the log A factor. Note that this is exactly the reason that the 
shattering threshold in the result of Barenboim et al. [BEPSvS] is 0(log^ A) rounds. 


3 The New Algorithm and Its Local Complexity 


Here we present a very simple and clean algorithm that guarantees for each node v that after 
0(logA + logl/e) rounds, with probability at least 1 — e, node v has terminated and it knows 
whether it is in the (computed) MIS or it has a neighbor in the (computed) MIS. 


The Intuition: Recall that the difficulty in locally analyzing Luby’s algorithm was the fact that 
the degree-dropping progresses of a node can be delayed by those of its neighbors, which in turn can 
be delayed by their own neighbors, and so on (up to log A hops). To bypass this issue, the algorithm 
presented here tries to completely disentangling the “progress” of node v from that of nodes that are 
far away, say those at distance above 3. 

The intuitive base of the algorithm is as follows. There are two scenarios in which a node v has 
a good chance of being removed: either (1) v is trying to join MIS and it does not have too many 
competing neighbors, in which case v has a shot at joining MIS, or (2) a large enough number of 
neighbors of v are trying to join MIS and each of them does not have too much competition, in which 
case it is likely that one of these neighbors of v joins the MIS and thus v gets removed. These two 
cases also depend only on v’s 2-neighborhood. Our key idea is to create an essentially deterministic 
dynamic which has these two scenarios as its (more) stable points and makes each node v spend a 
significant amount of time in these two scenarios, unless it has been removed already. 


The Algorithm: In each round t, each node v has a desire-level pt{v) for joining MIS, which 
initially is set to po{v) = 1/2. We call the total sum of the desire-levels of neighbors of v it’s 
effective-degree dt{v), i.e., dt{v) = The desire-levels change over time as follows: 

. fpt(?^)/2’ iidt{v)>2 

\min{2pt(u), 1/2}, if dt{v) < 2. 

The desire-levels are used as follows: In each round, node v gets marked with probability pt{v) 
and if no neighbor of v is marked, v joins the MIS and gets removed along with its neighbors'. 
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Again, each round of the algorithm can be implemented in 2 communication rounds on G, one 
for exchanging the desire-levels and the marks, and the other for informing neighbors of newly joined 
MIS nodes. Ignoring this 2 factor, in the sequel, each round means a round of the algorithm. 

The Analysis: The algorithm is clearly correct meaning that the set of nodes that join the MIS is 
indeed an independent set and the algorithm terminates at a node only if the node is either in MIS 
or adjacent to a node in MIS. We next argue that each node v is likely to terminate quickly. 

Theorem 3.1. For each node v, the probability that v has not made its decision within the first 
/3(logdeg -|- log 1 /e) rounds, for a large enough constant fi and where deg denotes v’s degree at the 
start of the algorithm, is at most e. Furthermore, this holds even if the outcome of the coin tosses 
outside {v) are determined adversarially. 

Let us say that a node u is low-degree if dfiu) < 2, and high-degree otherwise. Considering the 
intuition discussed above, we define two types of golden rounds for a node v: ( 1 ) rounds in which 
dt{v) < 2 and pt{v) = 1 / 2 , ( 2 ) rounds in which dv{t) > 1 and at least dt{v)/10 of it is contributed 
by low-degree neighbors. These are called golden rounds because, as we will see, in the first type, v 
has a constant chance of joining MIS and in the second type there is a constant chance that one of 
those low-degree neighbors of v joins the MIS and thus v gets removed. For the sake of analysis, let 
us imagine that node v keeps track of the number of golden rounds of each type it has been in. 

Lemma 3.2. By the end of round /3(log deg -|- log 1/e), either v has joined, or has a neighbor in, the 
(computed) MIS, or at least one of its golden round counts reached 100(logdeg -|- log 1/e). 

Proof. We focus only on the first /3(logdeg -|- log 1/e) rounds. Let gi and g 2 respectively be the 
number of golden rounds of types 1 and 2 for v, during this period. We assume that by the end of 
round /3(log deg -|- log 1 /e), node v is not removed and <71 < 100 (log deg -|- log 1 /e), and we conclude 
that, then it must have been the case that 52 > 100(logdeg -|- log 1/e). 

Let h be the number of rounds during which dfiv) > 2. Notice that the changes in pt{v) are 
governed by the condition dfiv) > 2 and the rounds with dfiv) > 2 are exactly the ones in which 
Pt{v) decreases by a 2 factor. Since the number of 2 factor increases in pt{v) can be at most equal 
to the number of 2 factor decreases in it, we get that there are at least /3(logdeg -|- log 1/e) — 2h 
rounds in which pt{v) = 1/2. Now out of these rounds, at most h of them can be when dfiv) > 2. 
Hence, gi > /3(logdeg -|- log 1/e) — 3h. As we have assumed gi < 100(logdeg -|- log 1/e), we get that 
/3(logdeg -|- log 1/e) — 3h< 100(log deg -|- log 1/e). Since (5 > 1300, we get h > 400(log deg -|- log 1/e). 

Let us consider the changes in the effective-degree dfiv) of v over time. If dfiv) > 1 and this is 
not a golden round of type- 2 , then we have 

dt+i{v) < 2^dfit) -h < \dt{v). 

There are 52 golden rounds of type-2. Except for these, whenever dfiv) > 1, the effective-degree 
dt{v) shrinks by at least a 2/3 factor. In those exceptions, it increases by at most a 2 factor. Each 
of these exception rounds cancels the effect of at most 2 shrinkage rounds, as (2/3)^ x 2 < 1. Thus, 
ignoring the total of at most 3g2 rounds lost due to type-2 golden rounds and their cancellation 


^There is a version of Luby’s algorithm which also uses a similar marking process. However, at each round, letting 
degiv) denote the number of the neighbors of v remaining at that time, Luby’s sets the marking probability of each 
node V to be which by the way is the same as the probability of v being a local minima in the variant described 

in Section 2. Notice that this is a very strict fixing of the marking probability, whereas in our algorithm, we change 
the probability dynamically/flexibly over time, trying to push towards the two desirable scenarios mentioned in the 
intuition, and in fact, this simple dynamic is the key ingredient of the new algorithm. 
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effects, every other round with dt{v) > 2 pushes the effective-degree down by a 2/3 factor®. This 
cannot (continue to) happen more than log 3 / 2 deg often as that would lead the effective degree to 
exit the dt{v) > 2 region. Hence, the number of rounds in which dt{v) > 2 is at most log 3/2 deg-|- 35 ( 2 - 
That is, h < log 3/2 deg-|- 35 f 2 - Since h > 400(log deg-|-log 1/e), we get (72 > 100(log deg-|-log 1/e). □ 

Lemma 3.3. In each type-1 golden round, with probability at least 1/100, v joins the MIS. Moreover, 
in each type-2 golden round, with probability at least 1/100, a neighbor of v joins the MIS. Hence, 
the probability that v has not been removed (due to joining or having a neighbor in MIS) during the 
first /3(logdeg + log 1/e) rounds is at most e. These statements hold even if the coin tosses outside 
N+{v) are determined adversarially. 

Proof. In each type-1 golden round, node v gets marked with probability 1/2. The probability that 
no neighbor of v is marked is nueAf('u)(^ — Pt{u)) > 4 “ ^ 4-2 _ Hence, 

V joins the MIS with probability at least 1/32 > 1/100. 

Now consider a type-2 golden round. Suppose we walk over the set L of low-degree neighbors 
of v one by one and expose their randomness until we reach a node that is marked. We will find a 
marked node with probability at least 

1 - ]J (1 -p„(t)) > 1 - > 1 - > 1 - > 0.08. 

u&L 

When we reach the first low-degree neighbor u that is marked, the probability that no neighbor 
of u gets marked is at least n«)eAr(n)(l ~ Pt{w)) > > ^-dtiu) y 1 /I 6 . Hence, with 

probability at least 0.08/16 = 1/100, one of the neighbors of v joins the MIS. 

We now know that in each golden round, v gets removed with probability at least 1/100, due to 
joining MIS or having a neighbor join the MIS. Thus, using Theorem 3.2, we get that the probability 
that v does not get removed is at most (1 — l/lOOj^'^^dogdeg+iogi/e) ^ e/deg < e. □ 

4 Improved Global Complexity 

In this section, we explain how combining the algorithm of the previous section with some known 
techniques leads to a randomized MIS algorithm with a high probability global complexity of 
0(log A) + 20 (Viogiogrx) rounds. 

As explained in Section 1 . 2 , the starting point is to run the algorithm of the previous section for 
0(log A) rounds. Thanks to the local complexity of this base algorithm, as we will show, we reach 
the shattering threshold after O(logA) rounds. The 2-hops randomness locality of Theorem 3.1, the 
fact that it only relies on the randomness bits within 2 -hops neighborhood, plays a vital role in 
establishing this shattering phenomena. The precise statement of the shattering property achieved 
is given in Theorem 4.2, but we hrst need to establish a helping lemma: 

Lemma 4.1. Let c > 0 be an arbitrary constant. For any 5-independent set of nodes S—that is, a 
set in which the pairwise distances are at least 5 —the probability that all nodes of S remain undecided 
after 0(clog A) rounds of the MIS algorithm of the previous section is at most 


^Notice the switch to dt(v) > 2, instead of dt{v) > 1. We need to allow a small slack here, as done by switching 
to threshold dt{v) > 2, in order to avoid the possible zigzag behaviors on the bonndary. This is because, the above 
argument does not bound the number of 2-factor increases in dt{v) that start when dt{v) £ (1/2,1) but these would 
lead dt{v) to go above 1. This can continue to happen even for an unlimited time if dt{v) keeps zigzagging around 1 
(unless we give further arguments of the same flavor showing that this is not possible). However, for dt{v) to go/stay 
above 2, it takes increases that start when dt{v) > 1, and the number of these is upper bounded to g 2 . 




Proof. We walk over the nodes of S one by one: when considering node v G S, we know from that 
Theorem 3.1 that the probability that v stays undecided after 0(clogA) rounds is at most A“'^, 
and more importantly, this only relies on the coin tosses within distance 2 of v. Because of the 
5-independence of set S, the coin tosses we rely on for different nodes of S are non-overlapping and 
hence, the probability that the whole set S stays undecided is at most □ 

From this lemma, we can get the following shattering guarantee. Since the proof is similar to that 
of [BEPSv3, Lemma 3.3], or those of [Bec91, Main Lemma], [ARVX12, Lemma 4.6], and [LRY15, 
Theorem 3], we only provide a brief sketch: 

Lemma 4.2. Let c be a large enough constant and B be the set of nodes remaining undecided after 
0(clog A) rounds of the MIS algorithm of the previous section on a graph G. Then, with probability 
at least 1 — 1 /n'^, we have the following two properties: 

(PI) There is no (G^ )-independent (G® )-connected subset S C B s.t. j5j > log^n. Here 
denotes the graph where we put edges between each two nodes with G-distance at most x. 

(P2) All connected components of G\B], that is the subgraph of G induced by nodes in B, have each 
at most 0 (log^ n • A^) nodes. 

Proof Sketch. Let H = G^ \ G^ , i.e., the result of removing G‘^~ edges from G®“. For (PI), note 
that the existence of any such set S would mean H[B] contains a (log^ n)-node tree subgraph. There 
are at most 4^°®^ different (log^ n)-node tree topologies and for each of them, less than nA*°®^ 
ways to embed it in H. For each of these trees, by Theorem 4.1, the probability that all of its 
nodes stay is at most By a union bound over all trees, we conclude that with probability 

1 — n( 4 A)^°s^ > 1 — l/n'^, no such such set S exists. For (P2), note that if G[B] has a 

component with more than 0(log^n • A"^) nodes, then we can find a set S violating (PI): greedily 
add nodes to the candidate S one-by-one, and each time discard all nodes within 4-hops of the newly 
added node, which are at most O(A^) many. □ 

From property (P2) of Theorem 4.2, it follows that running the deterministic MIS algorithm 
of Panconesi and Srinivasan [PS92], which works in rounds in graphs of size n', in each 

of the remaining components finishes our MIS problem in log rounds. However, the 

appearance of the log A in the exponent is undesirable, as we seek a complexity of O(logA) -|- 
20 (v 4 ogiogn^_ Xo remedy this problem, we use an idea similar to [BEPSv3, Section 3.2], which tries 
to leverage the (PI) property. 

In a very rough sense, the (PI) property of Theorem 4.2 tells us that if we “contract nodes that 
are closer than 5-hops” (this is to be made precise), the left over components would have size at 
most log^ n, which would thus avoid the undesirable log A term in the exponent. We will see that, 
while running the deterministic MIS algorithm, will be able to expand back these contractions and 
solve their local problems. We next formalize this intuition. 

The finish-off algorithm is as follows: We consider each connected component G of the remaining 
nodes separately; the algorithm runs in parallel for all the components. Eirst compute a (5, /i)-ruling 
set Rc in each connected component G of the set B of the remaining nodes, for an h = 0(loglogn). 
Recall that a (5, /i)-ruling set Rc means each two nodes of Rc have distance at least 5 while for 
each node in G, there is at least one node in Rc within its /i-hops. This (5, /i)-ruling set Rc can be 
computed in O(log log n) rounds using the algorithm® of Schneider, Elkin and Wattenhofer [SEW13]. 

®This is different than what Barenboim et al. did. They could afford to use the more standard ruling set algorithm, 
particularly computing a (5, 32 log A -|- 0(l))-ruling set for their purposes, because the fact that this 32 log A ends up 
multiplying the complexity of their finish-off phase did not change (the asymptotics of) their overall complexity. 
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See also [BEPSvS, Table 4], Form clusters around i^c-nodes by letting each node u G (7 join the 
cluster of the nearest iic'-node, breaking ties arbitrarily by IDs. Then, contract each cluster to a 
new node. Thus, we get a new graph G'^j on these new nodes, where in reality, each of these new 
nodes has radius h = O(loglogn) and thus, a communication round on G'(j can be simulated by 
0{h) communication rounds on G. 

From (PI) of Theorem 4.1, we can infer that G'q has at most log^n nodes, w.h.p., as follows: 
even though Rc might be disconnected in by greedily adding more nodes of G to it, one by 
one, we can make it connected in G‘^~ but still keep it 5-independent. We note that this is done only 
for the analysis. See also [BFPSv3, Page 19, Steps 3 and 4] for a more precise description. Since 
by (PI) of Theorem 4.1, the end result should have size at most log^n, with high probability, we 
conclude G'(j has at most log^ n nodes, with high probability. 

We can now compute an MIS of G, via almost the standard deterministic way of using network 
decompositions. We run the network decomposition algorithm of Panconesi and Srinivasan [PS92] on 
G'(j. This takes rounds and gives (7(^-clusters of radius at most logA*^)^ colored 

with 2 ®(\/^°siogA colors such that adjacent clusters do not have the same color. We will walk over 
the colors one by one and compute the MIS of the clusters of that color, given the solutions of the 
previous colors. Each time, we can (mentally) expand each of these G'^^ clusters to all the C-nodes 

of the related cluster, which means these C-clusters have radius at most log log n • "G ^ 

While solving the problem of color-j clusters, we make a node in each of these clusters gather the 
whole topology of its cluster and also the adjacent MIS nodes of the previous colors. Then, this 
cluster-center solves the MIS problem locally, and reports it back. Since each cluster has radius 
log log n • logA ^ takes log log n • logA rounds per color. Thus, over all the colors, 

the complexity becomes 2 ®^V^°siogA . Jog log n ■ 2 ®(V^°siogA = 2‘^(7iog logn) rounds. Including the 
O(loglogn) ruling-set computation rounds and the O(logA) pre-shattering rounds, this gives the 
promised global complexity of 0(log A) -|- 2 ‘^( 7 iogiogn)^ hence proving Theorem 1.2. 


5 Concluding Remarks 

This paper presented an extremely simple randomized distributed MIS algorithm, which exhibits 
many interesting local characteristics, including a local complexity guarantee of each node v ter¬ 
minating in 0(logdeg(u) -|- log 1/e) rounds, with probability at least 1 — e. We also showed that 
combined with known techniques, this leads to an improved high probability global complexity of 
0(log A) -|-2^(7iogiogn) rounds, and several other important implications, as described in Section 1.3. 

For open questions, the gap between the upper and lower bounds, which shows up when logA = 
u:{y/\og n), is perhaps the most interesting. We saw in (C2) of Section 1.3 that if the lower-bound 
is the one that should be improved, we need to go away from “tree local-views” topologies. An¬ 
other longstanding open problem is to hnd a poly(logn) deterministic distributed MIS algorithm. 
Combined with the results of this paper, that can potentially get us to an 0(log A) -|-poly(log log n) 
randomized algorithm. 

Acknowledgment: I thank Eli Gafni, Bernhard Haeupler, Stephan Holzer, Fabian Kuhn, Nancy 
Lynch, and Seth Pettie for valuable discussions. I am also grateful to Fabian Kuhn and Nancy 
Lynch for carefully reading the paper and many helpful comments. The point (C2) in Section 1.3 
was brought to my attention by Fabian Kuhn. The idea of highlighting the local complexity is rooted 
in conversations with Stephan Holzer and Nancy Lynch, and also in Eli Gafni’s serious insistence^^ 
that the (true) complexity of a local problem should not depend on n. 

^°This is a paraphrased version of his comment during a lectnre on Linial’s P(log* n) lower bound, in the Fall 2014 
Distributed Graph Algorithms (DGA) course at MIT. 
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A Simplified Global Analysis of Luby’s, due to Yves et al. 

We here explain (a slightly paraphrased version of) the clever approach of Yves et al. [YRSDZIO] 
for bounding Luby’s global time complexity: 

Lemma A.l. Let G[Vt] be the graph induced by nodes that are alive in round t, and let mt denote 
the number of edges of G\Vt]. For each round t, we have E[mt+i] < where the expectation is on 
the randomness of round t. 

Proof. Gonsider an edge e = {u,v) that is alive at the start of a round t, i.e., that is in ^[V)]- Note 
that edge e will not be in G[V)+i], in which case we say e died, if in the random numbers drawn in 
round t, there is a node w that is adjacent to u or u (or both) and w has the strict local minima of 
its own neighborhood. In this case, we say node w killed edge e. 
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Note that the probability that w kills e is , where d{w) denotes the degree of w. The 

difficulty comes when we want to compute the probability that there exists a w that kills e. This 
is mainly because, the events of different w killing e are not disjoint, and hence we cannot easily 
sum over them. Fortunately, there is a simple and elegant change in the definition, due to Yves et 
al. [YRSDZIO], which saves us from tedious calculations: 

Without loss of generality, suppose that w is adjacent to v. We say w strongly kills e from the side 
of node v —and use notation tc A e to denote it—if w has the (strictly) minimum random number 
in T^ UT^. Note that this is a stronger requirement and thanks to this definition, at most one node 
w can strongly kill e from the side of v. Thus, in a sense, we now have events that are disjoint which 
means the probability that any of them happens is the summation of the probabilities of each of 
them happening. The only catch is, we might double count an edge dying, because it gets (strongly) 
killed from both endpoints, but that is easy to handle; we just lose a 2-factor. In the following, with 
a slight abuse of notation, by E we mean the alive edges, i.e., those of G[V)], and by r(u), we mean 
the neighbors of v in G[V)]. We have 


K[Number of edges that die] > 


> 


> 


+ 


Pr[e gets strongly killed^ 

e=(v,w)£E 

E ( E Pr[t(; e] + ^ Pr[ri;' A ej'j /2 

e={v,w)£E hPEr(i;) 


( .d{w) + d{v) d{w)+d{u) 

e=(v^w)^E 


E E E 


+ E 


. X . . d{w) + div) ^ diw) + diu) 

vGV u^T{v}, e={u,v) ^w^T{v) w'&r{u) 
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E E 


d{w') + d{u) 

V) 
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d{v) d{w) 


^ ^ , d{w) + d{v) ' ^ d{w') -I- d{u) 

v^V tt;Er(rj) u£V w'£T(u) 


+ E E 


d(w) -|- d(v) ^ ^ diw) + div) 

v&Vw&T-{v) ^ ’ v&Vw&r{v) ^ ^ ^ ^ 


/2 


^ l)/2 = mt/2. 

vGV wGT{v} 


□ 


It follows from Lemma A.l that the expected number of the edges that are alive after 41ogn 
rounds is at most Therefore, using Markov’s inequality, we conclude that the probability 

that there is at least 1 edge that is left alive is at most Hence, with probability at least 1 — 
all nodes have terminated by the end of round 4 log n. 
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